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Abstract 

Temperature dependent measurements of ®^Fe Mossbauer spectra on CaFe 2 As 2 single crystals 
in the tetragonal and collapsed tetragonal phases are reported. Clear features in the temperature 
dependencies of the isomer shift, relative spectra area and quadrupole splitting are observed at the 
transition from the tetragonal to the collapsed tetragonal phase. From the temperature dependent 
isomer shift and spectral area data, an average stiffening of the phonon modes in the collapsed 
tetragonal phase is inferred. The quadrupole splitting increases by ~ 25% on cooling from room 
temperature to ~ 100 K in the tetragonal phase and is only weakly temperature dependent at 
low temperatures in the collapsed tetragonal phase, in agreement with the anisotropic thermal 
expansion in this material. In order to gain microscopic insight about these measurements we 
perform ah initio density functional theory calculations of the electric field gradient and the electron 
density of CaFe 2 As 2 in both phases. By comparing the experimental data with the calculations 
we are able to fully characterize the crystal structure of the samples in the collapsed-tetragonal 
phase through determination of the As z-coordinate. Based on the obtained temperature dependent 
structural data we are able to propose charge saturation of the Fe - As bond region as the mechanism 
behind the stabilization of the collapsed-tetragonal phase at ambient pressure. 
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I. INTRODUCTION 


CaFe 2 As 2 is conceivably one of the most interesting materials in the family of Fe-As 
based superconductors with a very fragile ground state and extreme sensitivity to the ex¬ 
ternal pressure and/or strain. At ambient pressure, a clear, sharp, hrst order transition 
at ~ 170 K from high temperature tetragonal - paramagnetic phase to low temperature 
orthorhombic - antiferromagnetic phase was observed in thermodynamic, transport and mi¬ 
croscopic measurements, in single crystals grown out of Sn flux.^^ Under pressure, this 
transition is initially suppressed at a rate of 10-15 K/kbai^^^^^. As the pressure increases, 
at about 3 kbar, a paramagnetic - collapsed tetragonal (cT) phase terminates the lower 
pressure orthorhombic - antiferromagnetic phase near 100 K. The cT transition tempera¬ 
ture increases under pressure and reaches 300 K at 15 kbar.l^ The uniaxial pressure 
derivatives of the structural/magnetic transition are large and of different signs.^The pres¬ 
sure - temperature phase diagram for CaFe 2 As 2 is very sensitive to non-hydrostaticity of the 
applied pressure 

This extreme sensitivity of pure, as well as substituted CaFe 2 As 2 single crystals to external 
stresses and strains yields an unexpected result: by a judicious choice of the annealing 
temperature and time (e.g. by control of internal stress helds via nanoscale precipitates) we 
can tune the structural/magnetic transition temperature in FeAs grown CaFe 2 As 2 crystals in 
a manner similar to application of pressure: suppressing the antiferromagnetic/orthorhombic 
transition by over 70 K and furthermore, we can obtain the crystals with the cT ground 
state at ambient pressure.^ By adding transition metal substitution as an extra parameter, 
we can tune the ground state of the Ca(Fei_xTa,) 2 As 2 (T = Co, Ni, Rh) crystals to be 
orthorhombic - antiferromagnetic, tetragonal - superconducting, tetragonal - paramagnetic, 
and collapsed tetragonal - paramagnetic, all at ambient pressurePressure measurements 
on Co-substituted CaFe 2 As 2 further illustrate the similar effects of applied pressure and post 
- growth annealing and quenching of sampled^ and also underscore the system’s pressure 
sensitivity with dT(^AF m/ ortho) / dP ~ —110 K/kbar and dTc/dP ~ —6 K/kbar. 

The aforementioned tunability of CaFe 2 As 2 created the opportunity for studies of the cT 
phase at ambient pressure (e.g. in FeAs grown CaFe 2 As 2 single crystals quenched from 960° 
C), using experimental techniques that are either difficult or impossible to combine with 
hydrostatic pressure.^^^® 
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^'^Fe Mossbauer spectroscopy has been widely used to study Fe-As based superconductors 
and related materials,including Sn-grown CaFe 2 As 2 single crystal^^^^^® and powdered 
single crystald^ primarily for probing magnetism (although it should be noted that grinding 
of CaFe 2 As 2 affects the macroscopic and microscopic magnetic propertie^^. In this work 
we will use this local probe technique to perform a detailed temperature-dependent study 
of the cT transition (and the cT phase) in single crystals of undoped, FeAs-grown 960° C 
quenchecP^, CaFe 2 As 2 . In these CaFe 2 As 2 samples both the high temperature tetragonal 
and low temperature cT phases are paramagnetic, with no magnetic order. We will compare 
the results with the available literature data on temperature and pressure induced transition 
to the cT phase. Some initial ^^Fe Mossbauer spectroscopy data on very similar FeAs-grown 
CaFe 2 As 2 crystals with cT ground state were published in Ref. [I5] and the pressure induced 
transition to the cT phase at room temperature in Sn- grown sample was recently studied 
by Mossbauer spectroscopy in Ref. [SH 

The isomer shift and the quadrupole splitting observed in Mossbauer spectroscopy can 
be directly related to the electron density at the nucleus and the electric field gradient 
respectively, both of which can be calculated in the framework of ab-initio Density Functional 
Theory. We will employ these calculations in order to provide a theoretical interpretation 
of the measurements, as well as to obtain a full structural characterization of the 960° C 
quenched samples in the collapsed-tetragonal phase by determining the As 2 ;-coordinates. 
We will also analyze the temperature dependent structural data in order to understand the 
microscopic mechanism responsible for the stabilization of the collapsed-tetragonal phase at 
ambient pressure. 


II. EXPERIMENTAL DETAILS AND COMPUTATIONAL METHODS 

CaFe 2 As 2 single crystals were grown from ternary melts rich in FeAs, following the pro¬ 
cedure described in Ref. US The excess liquid was decanted at 960° C, essentially quenching 
the samples from 960° C to room temperature. The cT transition temperature of the sample 
was ~ 90 K as conhrmed by dc magnetization measurements (see Fig. In the rest of 
the paper we will use CaFe 2 As 2 to refer to these FeAs grown, CaFe 2 As 2 samples quenched 
from 960° C to room temperature. 

Mossbauer spectroscopy measurements were performed using a SEE Co. conventional. 
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constant acceleration type spectrometer in transmission geometry with a ®^Co(Rh) sonrce 
kept at room temperatnre. The absorber was prepared as a mosaic of 0.04 - 0.1 mm thick 
single crystals held between two paper disks by a small amonnt of a pressnre sensitive 
adhesive. The gaps between the individnal crystals were kept as small as possible. The 
mosaic had the c - axis perpendicnlar to the disks (with the accnracy of or better than 
±10°) and arbitrary in-plane orientation. Since CaFe 2 As 2 and related materials are known 
to be very sensitive to the stresses and strains, inclnding those cansed by sample monnting,^ 
we performed dc magnetization measurements on a subset of the mosaic with the same 
mounting as for the Mossbauer measurements (Fig. [^. The cT transition is still sharp, 
without apparent shift of the transition temperature. 

The constrains on the thickness of the CaFe 2 As 2 crystals that can be cleaved without 
signhcant mechanical damage result in the absorber having 5-20 mg natural Fe per cm^ 
that corresponds to a ’’thick” absorber limit.l^The consequences of using ’’thick” absorber 
include thickness broadening of the absorption lines, increased error bars in absolute deter¬ 
mination of the quadrupole splitting and different degrees of saturation of the absorption 
lines of different intensity. The intrinsic linewidth found in inorganic compounds is often 
in the range of 0.26 - 0.28 mm/s,l^so the broadening observed in our measurements is not 
extremely severe. The other effects of the absorber thickness manifest as ’’smooth” contri¬ 
bution to the temperature dependent data and cannot cause the sharp features observed in 
the temperature dependences of the hyperhne parameters at the cT transition. 

The absorber holder consists of two nested white Delrin cups. The absorber holder was 
locked into thermal contact with a copper block with a temperature sensor and a heater, 
and was aligned with the 7 - source and detector. The c-axis of the crystals in the mosaic 
was parallel to the Mossbauer 7 beam. The absorber was cooled to a desired temperature 
using a Janis model SHI-850-5 closed cycle refrigerator (with vibration damping). The 
measurements were performed on warming, after initially cooling the absorber down to 
a base temperature of 4.6 K. The driver velocity was calibrated using an a-Fe foil and all 
isomer shifts (IS) are quoted relative to the a-Fe foil at room temperature. All the Mossbauer 
spectra were htted by the commercial software package MossWinn.l^ 

Ab initio Density Functional Theory (DFT) calculations were performed with the scalar 
relativistic Full Potential Local Orbital method (FPLO)l^. The electric held gradients were 
obtained according to the technique described in Ref. |39| and implemented in the FPLO 
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code. The Per dew-Burke- ErzenhoP^ exchange-correlation potential was employed for all 
calculations, and we considered 20 x 20 x 20 k-point meshes, with charge density converged 
within 10“® a.u.“^. In order to calculate the electron density at the location of atomic nuclei, 
we assumed the model of a hnite-size nucleus with constant charge density. Consistency 
of the calculated electric held gradients and electron densities was verihed against values 
obtained by the Full-Potential Linearized Augmented Plane-Wave method implemented in 
the WIEN 2 K codtP. 

III. EXPERIMENTAL RESULTS AND DISCUSSION 

The ®^Fe Mossbauer spectra of CaFe 2 As 2 at selected temperatures between the room 
temperature and the base temperature of 4.6 K are shown in Fig. At each temperature 
the spectrum is a quadrupole split doublet with the higher velocity line having almost twice 
the intensity of the lower velocity line. No secondary phase signal was observed in the 
spectra. Visually, the temperature induced changes, even near T^t ~ 90 K, are subtle. 

The temperature dependent linewidth of the spectra is shown in Fig. a). Overall 
the linewidth increases by a few percent on cooling from room temperature to the base 
temperature. As has been discussed in Ref. [151 fhe point symmetry and the location of 
the Fe site in the ThCr 2 Si 2 structure constrains the principal axis of the local electric held 
gradient tensor to the c - crystalline axis, as a result, a doublet lines intensity ratio of 3 : 1 is 
expected for the mosaic with the c - axis parallel to the 7 beam. Instead, the observed room 
temperature ratio is ~ 2.4 : 1 and it decreases to ~ 1.95 : 1 at the base temperature (Fig. 
(b)). The general behavior is consistent with the two temperature points reported earlier.^ 
Several possible reasons for the doublet lines intensity ratio being different from 3 : 1 were 
discussed in Ref. usi Most probably, an intensity ratio smaller than 3 : 1 and its further 
decrease on cooling may be due to the ’’thick” absorber conditions of the measurements 
(consistent with the large absorption observed (Fig. [^. We observe a small but distinct 
feature in the intensity ratio at the transition temperature around ~ 90 K. The change in 
the intensity ratio at the cT transition rehects an increase of the degree of saturation (caused 
by the thick absorber) of the stronger spectral line, thus pointing to stiffening of the lattice 
in the low temperature cT phase. 

In Fig. I^we show the measured isomer shift (IS) which increases upon cooling and has 
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a clear step-like change between 95 K and 85 K. The isomer shift includes contributions 
from both the chemical shift and the second-order Doppler shift. The latter is known to 
increase convexly upon decreasing temperature, due to gradual depopulation of the excited 
phonon states. However, it should be constant at low temperature, because of the quantum 
mechanical zero-point motion. The chemical shift should not depend on temperature. The 
main contribution to the temperature dependence of the isomer shift is therefore from the 
second-order Doppler shift, and is usually described by the Debye model®^ 


S{T) = 6{0) 


9kBT / T y x^dx 

io e--l’ 


( 1 ) 


where c is the velocity of light, M is the mass of the ®^Fe nucleus, and 5(0) is the temperature- 
independent part. For the isomer shift data in Fig. [^separate Debye hts, for the tetragonal 
phase (T > 95 K) and for the collapsed tetragonal phase (T < 85 K) were performed. The 
hts yield 0 _d = 239 ± 16 K for the high temperature tetragonal phase and 0 _d = 316 ± 7 
K for the low temperature cT phase. These results suggest that the lattice is stiffer in the 
cT phase. It has to be noted that the results of the ht in the low temperature cT phase 
should be taken with caution, since the ht is performed on a limited dataset with only 
minor temperature dependence. Still the conclusions are consistent with the rest of the data 
discussed in this work. 

The relative spectral area (Fig. also shows a distinct step - like feature in the region 
of the structural transition. Although we do not perform a Debye ht of the temperature 
dependent spectral area due to the possible ehects of the thick absorber, the step - like 
increase of spectral area (by ~ 6%) at the cT transition points to an increase of the Debye 
temperature, i.e. stihening of the lattice, in the cT phase. Somewhat higher, ~ 9%, increase 
in the relative spectral area is expected using 0 d values from hts of temperature-dependent 
isomer shift and the Debye model for the recoil-free fractionpSEH 


/ = exp 


- 3 ^^ 1 

kB&DMc^ 4 \0d/ 


i'Sd/T 


- 1 


( 2 ) 


(where / is the recoil-free fraction, and is the y-ray energy) leading to the conclusion 
that the value of Qb> in the cT phase may be slightly overestimated. 
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Stiffening of the lattice npon transition to a collapsed tetragonal phase (nnder pres¬ 
sure) was clearly observed in BaFe 2 AsJ^ via compressibility measurements, however similar 
set^SlMl of structural data for CaFe 2 As 2 appear to be less convincing. Nonetheless, a detailed 
study of phonon spectra in CaFe 2 AssP concluded that the phonons polarized in the ab 
plane are very similar in the tetragonal and cT phases, whereas transverse acoustic phonons 
propagating along the c direction stiffen very signihcantly in the cT phase. These results 
are consistent with our analysis of the changes in the Mossbauer spectra at the structural 
transition. 

The temperature dependent quadrupole splitting (QS) data is presented in Fig. The 
total change of the QS with temperature is ~ 25%. QS varies strongly in the tetragonal 
phase and is almost constant in the cT phase. The cT transition is clearly seen as a cusp in 
QS(T). Qualitatively, this behavior is consistent with the anisotropic thermal expansion of 
CaFe 2 AsJ 2 Il the quadrupole splitting is related to the c/a ratio, that changes considerably 
above the cT transition and signihcantly less in the cT phase. In a very general sense 
this behavior is consistent with the results from the Ref. |MI (with pressure as a control 
parameter): QS increases on approaching the cT transition, displays a cusp at the transition, 
and does not change in a signihcant manner in the cT phase. The origin of the cusp in the 
QS(T) data at the cT transition (possibly observed as a function of pressur^ as well) is 
however not clear and requires additional studies beyond the scope of the present work. 


IV. MISCROSCOPIC MODELLING 

To better understand the physics behind the experimental observations we rely on the 
fact that the isomer shift and the quadrupole splitting are directly related to the electronic 
charge densities and the electric held gradients, which can be calculated ab-initio in the 
framework of density functional theory (DFT). 

The isomer shift 5 can be expressecP^ as: 

5=^{Rl-Rl)[Pa{^)-Psm (3) 

O£o 

where Z is the nuclear charge, e is the elementary charge, eo is the permittivity in vacuum, 
Rg and Re are the nuclear ground and excited state radii and Ps(0) and Pa(0) are electron 
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densities at the source and absorption nuclei respectively. The only quantity that varies in 
the measurement is the electron density around the absorption nucleus Pa(0) so that we can 
simply write Eq. ([^ as a linear relation 


6 = a[pa(0) - b] 


( 4 ) 


Note that for ®^Fe the constant a is negativd^^®^. With this taken into account, part of the 
observed step-like increase in the isomer shift at the onset of the volume collapse (Fig. 
would correspond to the likewise drop in the electron density at the iron nucleus. This can 
be explained as an immediate consequence of Fe charge delocalization brought in by the 
volume collapse as a result of a strong, hrst order, structural phase transition. 

Furthermore, the quadrupole splitting A can be expressed a^^ 


A = 



( 5 ) 


where Q is the nuclear quadrupole moment, r] = \Vxx — Vyy\/\Vzz\ is the electric held gradient 
asymmetry, and \ Vxx\ A \^yy\ A iKzl are the eigenvalues of the electric held gradient tensor. 
Since the sample is tetragonal, the asymmetry p is zero and we obtain the following linear 
equation 


A = cVzz 


( 6 ) 


where c = eQ/2. 

Understanding the origin of the QS behavior (Fig. based on Eq. |^is quite involved since 
the electric held gradient Wz strongly depends on crystal structural details, in particular 
the details of the Fe-As bonding. 

The crystal structure of the tetragonal CaFe 2 As 2 is that of ThCr 2 SiJ^ and is described 
by the space group M/mmm. The Wyckoh positions corresponding to Ca and Fe atoms 
are fully hxed by the symmetry, while the Wyckoh position corresponding to the As atoms 
has only x and y coordinates hxed, leaving the z-coordinate as the sole internal structural 
degree of freedom. 



The structure of the 960° C quenched CaFe 2 As 2 sample was characterized in the 5-300 K 
temperature range by X-ray diffraction measurement^^. However, these measurements were 
able to determine only the unit cell parameters, and the As ^-coordinate was left unknown, 
rendering the crystal structure information incomplete. 

Full crystal structure information is the starting point for any DFT calculation. Since we 
don’t have here information of the As ^-coordinate but we do have valuable information from 
experiment on isomer shifts and quadrupole splittings, we are going to perform a “reverse 
engineering”, namely from the experimental data we will deduce the As z coordinates as a 
function of temperature. We would like to point out that the task of providing an accurate 
determination of the As ^-coordinates has additional importance through the fact that As 
atom positions are controlling parameters in the tetrahedral coordination geometry of the Fe 
atoms (see Fig. 04). and they are thus a very sensitive driver of physics of the iron-arsenide 

layeM 

Obtaining the As z-coordinate in a straightforward way, by structural optimization within 
DFT, may be rather inaccurate because it would require the accurate reproduction of the 
intricate internal mechanical conditions experienced by the samples used in the experiment. 
In this particular case, the sample strain, where Sn-grown single crystals are taken as a 
reference, is substantial enough to maintain the collapsed-tetragonal phase at the ambient 
pressure, and as a consequence, cannot be ignored in any theoretical description pretending 
to be an accurate portrayal of reality. 

To get around this problem, we chose a different strategy. We exploit the linear rela¬ 
tionship between the quadrupole splitting and the electric held gradient (Eq. ([^). The 
proportionality constant is the electron charge multiplied by the ®^Fe nuclear quadrupole 
moment, which is a well determined quantity, with values 0.15-0.17 barn®^®^. This allows 
us to compute the quadrupole splitting by calculating Vzz- We consider unit cell parameters 
experimentally determined at temperature calculate the quadrupole splitting Adft for 
a range of As 2 ;-coordinates and then perform a polynomial ht of the calculated data points 
to obtain a curve Adft(7’, z(T)), where z(T) denotes As ^-coordinate at temperature T. We 
then determine the ^-coordinates by solving numerically Adft(T, z(T)) = Aexp(T), where 
Aexp(T) is the quadrupole splitting experimentally determined at temperature T. The As 
^-coordinates obtained in this way are shown in Fig. along with the height of As atoms 
above the Fe plane. Although the quadrupole splitting shows non-monotonous behavior 
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in the 0-90 K region, reversing the decrease around 30 K, the As ^-coordinate, and corre¬ 
sponding height above the Fe plane display a monotonous increase and decrease respectively, 
indicating the well known complex relationship between the position of the As atom and 
the electronic propertied. 

Since the As ; 2 -coordinates now satisfy Eq. ([^ by construction, being able to simulta¬ 
neously satisfy Eq. Q provides an important consistency check. In Fig. we show the 
relationship between the measured isomer shift and the electron density calculated for the 
As 2 ;-coordinates obtained from the quadrupole splitting. The data is htted very well by a 
linear function and the slope is determined to be a = —0.45 ± 0.04 [a.u.^ mm/s] which is in 
good agreement with previously reported values for ^^Fe in various compoundd®^. 

It is important to note here that Figs. and [8] are restricted to the collapsed-tetragonal 
phase, since in the tetragonal phase the experimentally measured values of the quadrupole 
splittings cannot be reproduced by the calculations within a reasonable range of As 2 ;- 
coordinate values. We attribute this to the fact that antiferromagnetic fluctuations and 
electronic correlation effects cannot be accurately captured within DFT. In the collapsed 
tetragonal phase it has been shown that both correlation effectd^^^^ and antiferromagnetic 
fluctuation^® are suppressed resulting in a good agreement between our DFT calculations 
and the Mossbauer measurements. 

Focusing now on the tetragonal phase, it is well known that matching the experimentaly 
determined! and the DFT calculated equilibrium. As ^-coordinates in Fe-based supercon¬ 
ductors is, in general, problematic as has been pointed out by numerous publicationd^®^. 
However, spin polarized calculations within GGA have been shown to provide reasonable 
estimates, specially for the orthorhombic phases with long range magnetic ordei®. 

Such approach can be used here to mimic the effect of electron correlation and antifer¬ 
romagnetic fluctuations in the tetragonal phase in the framework of DFT. We have experi¬ 
mented with spin-polarized GGA calculations with and without on-site Goulomb repulsion. 
In this way we were able to extract more reasonable estimates for As ^^-coordinates. Here 
we present results with G=3 eV and Hund’s coupling J=1 eV. We assume Neel order in 
order to preserve the tetragonal symmetry and account for the antiferromagnetic character 
of fluctuations. In this way we obtain an estimate of the As ; 2 -coordinate in the tetragonal 
phase which should be taken only as an indication of the general temperature trend. An 
LDA-I-DMFT (Local Density Approximation plus Dynamical Mean Field Theory) calcula- 
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tion would be more appropriate but it is beyond the scope of the present work. 

Figure]^ a) shows the DFT predicted Fe-As bond lengths both in the collapsed tetragonal 
and the tetragonal phase. As the sample is cooled down the Fe-As bond length decreases, 
driven by thermal contraction. Interestingly, this trend reverses in the collapsed tetragonal 
phase where the Fe-As bond length increases, working against the thermal contraction of 
the sample. It is instructive to compare this behavior to the experimentally determined 
Fe-As bond lengths at ambient pressure and 0.63 GPa in Sn-grown samples subjected to 
hydrostatic pressur^. We would like to stress here that the data shown corresponds to 
available measurement^® for CaFe 2 As 2 performed under different conditions to those used in 
the present study and thus our ab initio DFT results are not meant to reproduce these data. 
However, as we will show below, comparison to Sn grown, low strain, samples at P = 0 and 
P = 0.63 GPa provides references to help dehne the expected range of behavior seen in our 
GaFe 2 As 2 sample. 

All measurements in our study are performed at ambient pressure and on samples that 
have a certain amount of strain relative to the Sn-grown samples from Ref. 0 We can 
thus consider the samples used in our study to be intrinsically under certain amount of 
pressure, even at zero external pressure. Keeping this analogy in mind we can proceed now 
with the structural comparison. We observe that our predicted bond lengths £t into the 
experimentally observed range of bond lengths and that, despite the fact that the predicted 
bond lengths in the tetragonal phase are just a rough estimate, they are consistent with 
the experimental data when we take into account the differences in experimental conditions. 
Namely, the hydrostatic pressure stabilizes the cT phase below 200 K at 0.63 GPa in Sn- 
grown samples, while in our 960° G quenched samples the cT phase is stabilized below 90 
K and at ambient pressure. Hence, the ambient-pressure 960° G quenched sample structure 
should roughly correspond to the structure of the Sn-grown sample at an intermediate 
pressure between ambient pressure and 0.63 GPa which is precisely what the Fe-As bond 
lengths in Fig. [^a) show. The exception is a 50 K interval around the T-cT transition which 
is a consequence of hysteresis effects associated with the volume collapse. Since the collapsed 
phase is much stiffer than the tetragonal phas^^®^®^ we expect reduction of the structural 
differences between the Sn-grown samples under pressure and the 960° G quenched samples, 
which is again conhrmed by our data. Similar consistency in bond length values as well as 
jump at Tct is observed in the behavior of interlayer As-As bonds, shown in Fig. [^b). 
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However, the most striking feature observed in Fig. a) is the reversal of Fe-As bond 
contraction in the cT phase and we believe that this is the key to understanding the mech¬ 
anism behind the stabilization of cT phase at ambient pressure. Since, the formation of the 
cT phase is precipitated by the formation of interlayer As-As bonds, we need to look closer 
into the interplay of Fe-As and As-As bonding dynamics (Fig. [^b)). 

Examination of the behavior of the V^z electric held gradient component for the iron 


and arsenic nuclei (Fig. 10) in the collapsed tetragonal phase in CaFe 2 As 2 shows that 
each component follows a different trend. While the electric held gradient decreases at the 
iron nucleus, it increases at the arsenic nucleus as a consequence of the bonding behavior. 
Relative weakening of the Fe-As bond implies that charge density in the bonding region 
is reduced, leading to the reduction of the electric held gradient on the iron site. This is 
rehected in the angular momentum resolved electric held gradient tensoi^^ which shows a 
major contributions coming from the I = 1 component, implicating reduced hybridization 
of iron 3d orbitals with arsenic 4p orbitals as a main culprit leading to the reduction of 
the electric held gradient. On the other hand, the increase of the electric held gradient at 
the arsenic site follows from the stronger interlayer As-As bonding as the sample is cooled. 
This is also rehected in the angular momentum resolved electric held gradient, with I = 1 
component being essentially the only contribution. Such bond dynamics indicates that some 
of the charge in the Fe-As bond region is transfered to the interlayer As-As bond region, 
contributing to its strength and helping to stabilize the collapsed tetragonal phase at ambient 
pressure. 

For further answers we examine in Fig. j^b) the predicted interlayer As-As bond lengths 
for the 960° C quenched sample at ambient pressure. We again compare with the ehect 
of the hydrostatic pressure and show the measured values for Sn-grown CaFe 2 As 2 samples 
at the pressure of 0.63 GPa. We observe that the As-As bonds appear to be relatively 
longer in the tetragonal phase of the 960° C quenched sample. Such bond conhguration 
can lead to charge saturation in the Fe-As bond region under thermal contraction, meaning 
that at a certain temperature it becomes energetically more favorable to transfer some of 
the charge into the relatively less charge dense interlayer As-As bond region, precipitating 
the formation of the interlayer As-As bond which in turn stabilizes the collapsed tetragonal 
phase. This process continues in the collapsed tetragonal phase, where the Fe-As bond 
keeps getting weaker while the interlayer As-As bond gets stronger until zero temperature is 
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reached, where the bond lengths revert to the ratio observed in the Sn-grown samples under 
hydrostatic pressure. 


V. SUMMARY 

Experimental clear features in the temperature dependencies of the isomer shift, relative 
spectral area and quadrupole splitting are observed at the transition from tetragonal to 
cT phase in CaFe 2 As 2 . From the analysis of the temperature dependent isomer shift and 
spectral area data, as well as the temperature evolution of the intensity ratio, an average 
stiffening of the phonon modes in the cT phase is inferred. Quadrupole splitting increases 
by 25% on cooling from room temperature to the cT transition, in agreement with a 
large anisotropy of thermal expansion in this temperature region. A clear unusual cusp was 
observed in the QS(T) data at the cT transition. The origin of it is not clear and requires 
additional studies. 

Additionally, we performed ab initio density functional theory calculations of the electric 
field gradient and the electron density of CaFe 2 As 2 in the tetragonal and collapsed tetragonal 
phase. Comparison of these calculations to the measurements allows for the determination of 
the As ; 2 -coordinates in the 960° C quenched sample, providing the full structural information 
under these conditions and allowing us to probe the microscopic origin of the experimental 
observations. 

We have found that the resulting Fe-As bond lengths are consistent with the data recorded 
for Sn-grown samples under hydrostatic pressure, as expected from the broad similarity of 
the observed structural effects and proper accounting of differences in experimental condi¬ 
tions. In the tetragonal phase, as the sample is cooled, the Fe-As bonds undergo thermal 
contraction, resulting in larger electric held gradients at the iron nuclei and leading to an 
increase of the quadrupole splitting. In the collapsed tetragonal phase we observe instead 
an anomalous behavior of the Fe-As bond, whereby it undergoes expansion as the sample is 
cooled, resulting in a decrease of the observed quadrupole splitting. We attribute this be¬ 
havior to charge redistribution from the Fe-As bond region into the interlayer As-As region, 
increasing the As-As bonding strength and resulting in stabilization of the collapsed tetrag¬ 
onal phase. In our proposed scenario, it is high temperature quenching at 960° C which 
stabilizes the relatively shorter Fe-As bonds, leading to charge saturation of the Fe-As bond 
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region as the sample contracts. 
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FIG. 1. Temperature dependent magnetization, M/H{T) and its temperature derivative, 
d{M/H)/dT, of a part of CaFe 2 As 2 crystals mosaic used for Mossbauer measurements in the 
vicinity of the cT transition. Note that the paper and adhesive of the mount give some small, 
weakly temperature dependent contribution to M/H{T). 
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FIG. 2. (Color online) Mossbauer spectra of CaFe 2 As 2 at selected temperatures. Symbols - 
data, lines - fits. T^t ~ 90 K, as shown in Fig. 
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FIG. 3. (Color online) CaFe 2 As 2 ®^Fe Mossbauer data: (a) Temperature dependence of the 
linewidth; (b) The intensity ratio of the doublet lines as a function of temperature. Circles - 
this work, rhombuses - Ref. ESI Dashed lines are linear fit to the data [panel a)] and guides to the 
eye [panel b)]. 
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FIG. 4. (Color online) Temperature dependent isomer shift of CaFe 2 As 2 . Dashed line - Debye fit 
for T > 95 K resulting in 0£) = 239 ± 16 K, dash-dotted line Debye fit for T < 85 K resulting in 
0£, = 316±7 K 
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FIG. 5. (Color online) Temperature dependent relative spectral area of CaFe 2 As 2 . Dashed lines 
are guides to the eye. 6% change is shown as a vertical bar. 
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FIG. 6. (Color online) Temperature dependent quadrupole splitting of CaFe 2 As 2 . Circles - this 
work, rhombuses - Ref. ESI 
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FIG. 7. (Color online) (a) Calculated temperature dependence of As z-coordinate in the collapsed- 
tetragonal phase. The inset shows temperature dependence of height of As atoms above the Fe 
plane, (b) Diagram showing geometry of the FeAs layer. 
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FIG. 8. (Color online) The linear relationship between the calculated charge density at the iron nu¬ 
cleus and the measured isomer shift in the collapsed-tetragonal phase. Note: implicit temperature 
variation is shown by temperature values for the end points. 
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FIG. 9. (Color online) (a) The ab-initio predicted Fe-As bond length in the 960° C quenched sample 
compared against the measured bond lengths on Sn-grown samples under hydrostatic pressur^. 
Dotted lines are guides for the eye. (b) Predicted interlayer As-As bond lengths in the 960° C 
quenched sample compared against the Sn-grown sampleJ^. Vertical line marks Tct at P = 0.63 
GPa. Please note that the comparison to experimental data is intented to estimate the properties 
of the “intrinsic” pressure in 960° C quenched sample (see the text for details). 
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FIG. 10. (Color online) The calculated 14^ component of the electric held gradient for the iron 
and arsenic sites normalized to their respective zero temperature values. 





